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Abstract 

The method of current extraction under linear increasing voltages (CELIV) allows for the si- 
multaneous determination of charge mobilities and charge densities directly in thin films as used 
in organic photovoltaic cells (OPV). In the past, it has been specifically applied to investigate the 
interrelation of microstructure and charge transport properties in such systems. Numerical and 
analytical calculations presented in this work show that the evaluation of CELIV transients with 
the commonly used analysis scheme is error prone once charge recombination and, possibly, field 
dependent charge mobilities are taken into account. The most important effects are an apparent 
time-dependence of charge mobilities and errors in the determined field dependencies. Our results 
implicate that reports on time-dependent mobility relaxation in OPV materials obtained by the 
CELIV technique should be carefully revisited and confirmed by other measurement methods. 
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I. INTRODUCTION 



The understandig of charge transport properties of organic semiconducting materials is 
one of the keys for further improvement of devices such as organic light-emitting diodes 
and organic photovoltaic cells. Owned to it's simplicity, straightforward data analysis and 
applicability to measurements on films of well below 100 nm thickness, the technique of cur- 
rent extraction by linear increasing voltages (CELIV, see Ref. jl!) has attracted considerable 
interest over the past few years. It's unique attractiveness stems from the opportunity to 
study charge transport directly in thin film systems as used in the actual devices. As such, it 
is the ideal technique to investigate structure-property relationships of charge transport and 
recombination e.g. in state-of-the-art donor /acceptor photovoltaic systems such as poly- 
mer/polymer and polymer/small molecule material blends, the morphology of which cannot 
be well reproduced in modified sample geometries. CELIV has been originally introduced to 
determine charge mobility and concentration in microcrystalline Si:H semiconductors and 
doped conjugated polymers under thermal equilibrium conditions,™ 1 ^ but was later applied to 
study the field-dependent mobility of non-equilibrium charge carriers created by the absorp- 
tion of a light pulse.--— Hereby photo-CELIV experiments have been used to quantify the 
decay of carrier density with time and to measure the time dependence of mobility. Time- 
dependent relaxation of charge mobilities as determined by the photo-CELIV technique has 
been reported for various state-of-the-art photovoltaic bulk hetero junction blends.--— -—^ It 
has usually been attributed to the relaxation of carrier energies in an extended density of 
transport states present in these disordered solids. 

This paper treats the analytic and numerical analysis of CELIV experiments under con- 
ditions diverging from the equilibrium assumptions of its original derivation. We first review 
the basic theory of equilibrium CELIV experiments and point out inaccuracies in the origi- 
nal derivations. Based on this, we treat the case of non-equilibrium photogenerated charge 
carriers with constant mobility by both numerical as well as analytical methods. We show 
that with the standard evaluation scheme applied under this condition, an apparently time- 
dependent charge mobility is deduced that is solely the effect of the used analysis. We 
provide an analytic analyis method to determine charge mobilities in a reliable way even 
under conditions of non-equilibrium charge carrier recombination. Finally, we treat the ap- 
plication of CELIV measurements to the situation of field-dependent carrier mobilities. A 
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modified analysis method is introduced and shown to result in significantly reduced errors 
for the derived field dependence. 



II. CHARGE MOBILITY DETERMINATION BY THE CELIV TECHNIQUE 

The experimental procedure and analytic evaluation for CELIV experiments has been 
introduced by Juska et al. in Ref. oJ, to which the reader is referred for further details. 
Consider the extraction of equilibrium charge carriers of density n and mobility fi in the 
electric field U/d, where d is the layer thickness and U(t) = U't is the applied voltage that 
rises linearly in time. Without restriction of the generality of our analysis, we assume that 
the internal electric field is zero at U = 0. This differs from the more general situation 
where electrodes of different work functions are used. Note that experimentally one should 
account for the built-in potential due to differences in the electrode work functions to ensure 
that the electric field within the organic layer is zero at U = 0. Assuming that the electrode 
area A is much larger than d, that one carrier type is much more mobile than the other 
(here: holes) and that the electrodes are non-injecting, the time dependent charge density 
is p(z,t) = —en for < z < l(t) and p(z,t) = for l(t) < z < d. Hereby, all holes have 
been depleted from the layer up to the time-dependent extraction depth l(t). The current 
density measured in the external circuit due to the extraction of charges at z — d is 
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(l(t)<d) 
else, 

assuming t ^> RC where R is the external circuit resistance and C = eA/d is the geometrical 
sample capacitance assuming a permittivity of e = e r eo- The extraction depth is the solution 
of 

dl(t) enpl 2 (t) _ fiU't 
dt + 2de ~ d ^' 
under the initial conditions 1(0) = and dl(0)/dt = 0. This is a first order nonlinear 
differential equation of Ricatti type that can be solved numerically parametric in the di- 
mensionsless voltage slope e 2 U' /2e 2 n 2 fid 2 . At some time t max) the current density ([T]) will 
peak at j(t max ), where the relative height Aj/j = (j(t mSuX ) — j(0))/j(0) can be expressed 
as a bijective function of the dimensionless parameter x — f i U't 2 nax /2d 2 . Figure [1] shows 
the result of our numerical calculation of x as a function of Aj/j(0). \ i s equal to 1/3 at 
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FIG. 1: Results of the calculation of x as function of Aj/j(0) from the numerical solution of 
equations (pQ) and (j2J) (symbols), compared to fit I and II as discussed in the text as well as the 
parametrisation used by Juska et al. in Ref. [2j. The inset shows a typical CELIV current calculated 
for d = 100 nm, U' = 2 x 10 4 V/s, fj, = 2 x 10 6 cm 2 /Vs, n = 10 22 m" 3 , A = 1 mm 2 and e r = 2.9, 
indicating the capacitive charging current j(0) and the maximum current j(^max)* 

Aj/j(0) = and decreases in a nonlinear fashion for Aj/j(0) > 0. Also shown in figured] 
are two different parametrizations of this curve in terms of \ — [3(1 + O.lSAj/j^O))]" 1 (fit 
I) and x = 0.329 exp[-0.180Aj/j(0)] + 0.005 exp[0.253Aj/j(0)] (fit II) which are good ap- 
proximations for Aj/j(0) < 1 or Aj/j(0) < 7 for fit I and fit II, respectively. Using these, 
the charge mobility can be calculated as ji = 2d 2 x/U't^ x , resulting in 
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when using fit II valid for Aj/j(0) < 7. This is in variance with the result published 
by Juska et al.— " and used in several articles,-* 1 ^ 3 - which corresponds to choosing x = 
[3(1 + 0.36Aj'/j(0))] _1 . In order to provide an independent test for the consistency of 
our analysis, we simulated the CELIV experiment with a numerical drift-diffusion solver 
program. Figure [2] shows the resulting current transients assuming a charge mobility of 
fi = 2 x 10~ 6 cm 2 /Vs, a layer thickness of d = 65 nm with e r = 3, a voltage slope of 
U' = 2 x 10 4 V/s and an initial charge carrier density of n = 4 x 10 22 m -3 . As expected, the 
numerically evaluated solution to equations (JTl) and (j2J) closely follows the simulation results 
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FIG. 2: Numerically simulated CELIV transients (solid lines), for simulation parameters see text. 
The results are shown both with and without charge diffusion according to the Einstein relation. 
The dashed line corresonds to the numerical solution of equations ([I]) and ([2]) and closely follows 
the simulation results excluding charge diffusion. 

when charge diffusion is suppressed. We determined Aj/j(0) = 0.781 and t max = 25.125 /is 
from the simulation data, from which the apparent mobilities /i = 1.96 x 10 6 cm 2 /Vs and 
fi = 1.95 x 10 6 cm 2 /Vs are calculated using equation ([3]) or equation (jl]), respectively. Using 
the Juska et al. result, we instead obtaine = 1.74 x 10~ 6 cm 2 /Vs, which proves that this 
approximation underestimates the mobility. We therefore suggest to calculate the mobility 
using equation (J3j) or (@J, depending on the magnitude of Aj/j(0). 

III. THE ROLE OF BIMOLECULAR CHARGE RECOMBINATION 

The gaussian disorder model of charge transport in organic semiconductors predicts that 
after photoexcitation, charge carriers will relax energetically towards their equilibrium en- 
ergy, with a concomitant mobility decrease. The understanding of this process is of consider- 
able importance for organic photovoltaic devices in order to further improve their efficiency. 
Recently, the CELIV method has been applied to the study of non-equilibrium charge car- 
riers photogenerated by the absorption of short laser pulses.- 1 ^ Various publications used 
this photo-CELIV technique to study both mobility and density relaxation of photogener- 
ated carriers .-^ii 2 - 1 ^ We feel that it is important to point out that the analytic equations 
used to evaluate CELIV experiments have been derived by assuming the presence of equi- 
librium carriers. Thus, bimolecular charge recombination has not been taken into account 
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FIG. 3: Comparison of numerically simulated CELIV current transients without (left) and with 
(right) bimolecular charge recombination according to the Langevin mechanism, for simulation 
parameters see text. The results are shown parametric in the initial charge density, the position of 
j (imax) is marked by a circle and vertical line. 

for the calculation of the time- dependent charge density p(z,t). However, assuming charge 
recombination according to the Langevin mechanism^ 4 -, the charge density at z > l(t) de- 
cays as n(t) = n(0)(l + t/r^) -1 , where r a = e/en(0)fi is the dielectric relaxation time. 
This renders the evaluation of CELIV experiments in terms of equation (j4]) inaccurate and 
motivates to study the effects of charge recombination by numerical simulations. Figure [3] 
shows the results of numerically simulated current transients assuming fi — 2 x 10 6 V/m, 
d = 100 nm, e r = 3 and U' — 2 x 10 4 V/s while varying the initial charge density n(0) 
between 10 22 m -3 and 10 23 m -3 . When bimolecular charge recombination is taken into ac- 
count in the simulations, a significantly reduced amount of charges is extracted, reducing 
Aj/j(0) and strongly shifting t max towards shorter times. Thus, the apparent mobilities 
calculated from transients affected by recombination using equation are expected to be 
higher compared to those calculated from recombination-free transients. We analyzed this 
in more detail for polymer blends of poly[2,5-dimethoxy-l,4-phenylenevinylene-2-methoxy- 
5-(2-ethylhexyloxy)-l,4-phenylenevinyiene] (M3EH-PPV) with poly[oxa-l,4-phenylene-l,2- 
(l-cyano)-ethylene-2,5-dioctyloxy-l,4-phenylene-l,2-(2-cyano)-ethylene-l,4-phenylene] (CN- 
ether-PPV), for details of these materials see Ref. ^. Samples were fabricated by spin- 
coating a 1:1 blend of these polymers from chlorobenzene solution onto precleaned and 
structured ITO substrates covered by a layer of PEDOT:PSS (Clevios AI4083 obtained 
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FIG. 4: CELIV current transients measured for a 55 nm thick M3EH-PPV:CN-ether-PPV blend 
layer at U' = 1.06 V/us for different delay times after photoexcitation ranging from near zero up 
to 500 /us. Connected circles indicate the determined j(t max ) points. The inset schematically shows 
the layer structure and illumination direction used for the experiment. 

from H.C Starck, Germany) and evaporating a 200 nm thick aluminum top electrode. De- 
vices were fabricated under protective nitrogen atmosphere and encapsulated by a cover 
glass and two-component expoxy resin prior to measurements under ambient conditions. 
Figure H] shows photo-CELIV current transients obtained for a device with a 55 nm thick 
polymer layer and an electrode area of A = 1 mm 2 at a voltage slope of U' = 1.06 V/ytis 
for various delay times td between photogeneration using a 20 ns long laser pulse of 355 nm 
wavelength and the beginning of charge extraction. As is obvious from the current tran- 
sients, the time t max of maximum extraction current strongly shifts to smaller values for 
short td- We modeled the impact of charge density on these measurement results by numer- 
ical simulation of the experiment using the model parameters d — 55 nm, U' = 1.06 V/jits, 
e = 3 and assuming a field- and time-independent mobility fi = 3.8 x 10 -6 cm 2 /Vs. The 
charge density n(0) at the beginning of charge extraction was varied between 10 21 m~ 3 and 
10 25 m -3 . Figure [5] shows the Aj/j(0) and apparent mobilities calculated from the sim- 
ulated CELIV transients using equation (j4j). The apparent mobility (solid symbols) rises 
with charge density for n(0) > 10 23 m -3 , corresponding to current maxima (open symbols) 
of Aj/j(0) > 0.1. It has been suggested^ that CELIV transients are most convenient 
to determine experimentally when Aj/j(0) ~ 1. Our simulations strongly discourage this 
choice whenever nonequilibrium charge carriers are investigated, since charge recombination 
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FIG. 5: Analysis of numerically simulated CELIV transients in terms of current maximum Aj/j(0) 
(open symbols) and apparent mobility (solid symbols) calculated from equation ([3]) for a range of 
charge densities n(0) present at the beginning of charge extraction (id = 0). The simulation model 
parameters correspond to the experimental situation of figure 01 

strongly distorts the transients in this regime. The simulation data can also be calculated 
by varying the delay time t d between photogeneration of charge carriers of density n p h to 
and charge extraction, whereby the charge density at the beginning of charge extraction is 
n(0) = n p hoto/(l + e/m photo t/e). Figure [6] compares the apparent mobilities calculated from 
simulation data using equation (J3J) with those obtained in the same way from figure HI Under 
the assumption of n p hoto = 4 x 10 23 m~ 3 , simulated results closely follow those obtained from 
the measurement over the whole range of delay times. For the shortest delay time, we esti- 
mate from the experimental current transient that a total amount of charges corresponding 
to a charge density of only 5 x 10 22 m -3 could be extracted during the CELIV experiments. 
For the simulated transients at a comparable delay time of = 0.38 /xs we calculated a very 
similar value of 6 x 10 22 m -3 . Thus, more than 80% of the initially generated charge pairs 
are lost due to bimolecular recombination during the delay time and the initial part of the 
CELIV transient. 

Since the effect of charge recombination during non-equilibrium CELIV experiments is 
of considerable experimental interest, we will elaborate further on an analytic treatment of 
this situation. In general, the charge density p will follow a complicated spatial and time 
dependence, since charge recombination will take place only in the region z > l(t), where 
both charge types are present. To keep our analysis sufficiently general, we consider the case 
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FIG. 6: Apparent mobility calculated from the experimental CELIV transients shown in figure H] 
using equation (symbols), compared with the results of numerical simulations assuming various 
photogenerated charge carrier densities n p hoto) a time- independent charge mobility /i = 3.8 x 
10 -6 cm 2 /Vs and bimolecular recombination according to the Langevin mechanism. Data is shown 
as a function of the delay time between photoexcitation and the beginning of the extraction 
voltage pulse, the approximate uncertainty of delay time for first data point is indicated by error 
bars. 

of a dielectric relaxation time given by r CT = e/en(0)/i/3, where ft is a recombination prefactor 
{(3 = 1 corresponds to Langevin recombination). Reduced bimolecular recombination with 
/3 < 1 has been shown to prevail in some polymer/small molecule donor /acceptor blend 
solar cell materials.— The spatial charge distribution for z <l{t) then depends on l{t) as 

-en 



P{z) 
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where /(/^^(z)) = z defines the inverse function Z( _1 )(z) of l(t). Despite this implicit 
expression for the charge density, analysis shows that the current transient can be obtained 
by the surprisingly simple expression 

1 - l/d dl(t) 
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The time dependence of l(t) is calculated from 

d 2 l{t) pU' 1 



l(t) dl(t) 
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dt 2 d d(3T a 1 + t/r a dt ' 
the solution of which converges to that of equation ([2]) for ft — > 0. We compared this solution 
to the results of numerical simulations and found good agreement when charge diffusion was 
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FIG. 7: Results of the calculation of x(fl)/x(fy as function of Aj/j(0) from the numerical solution 
of equations ([6]) and ((JJ), parametric in the bimolecular recombination prefactor /?. 

neglected in the simulations. Using the same techniques as for figure CD, we calculated x as 
a function of Aj/j(Q) parametric in the prefactor (3. Figure [7] plots the results relative to 
X as determined for (3 = 0. Given a specific recombination prefactor /3, these curves can be 
used to directly extract the recombination-corrected charge mobility from the measurement 
or to estimate the impact of recombination on CELIV results obtained by the more simple 
estimates provided by equations ([3]) and (HI). We additionally fitted x as a function of 
Aj/j(0) for the case (3 = 1 by the double exponential expression used above. This directly 
results in the expression 



which is valid at Aj/j(0) < 0.95 with a relative error of less than 3.5% and can be used 
to determine true charge mobilities from CELIV transients even under conditions of high 
charge densities, assuming that Langevin recombination prevails. 

IV. FIELD-DEPENDENT CHARGE MOBILITIES 

The charge mobility in organic semiconductors is usually considered to be both a field 
and density dependent quantity, where the electric field dependence has been experimentally 
found to mostly follow a Poole- Frenkel type law n{E) = /jLq exp(n\^E) . The work of Bassler et 
al. has shown that this can be understood in terms of transport sites having random energies 
according to a gaussian distribution, rendering /x and k temperature dependent.— CELIV 
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experiments provide a unique opportunity to determine charge mobilities in undoped organic 
semiconductor films of well below 100 nm thickness, but has the disadvantage of working 
under conditions of a non-constant electric field. Since the mobility is calculated from the 
maximum extraction current point, the determined values have usually been associated with 
the electric field E* = U't max /d (extraction field) present at the time of maximum extraction 
current although the validity of this approach has never been tested rigorously. As we 
have shown above, mobilities determined for non-equilibrium charge carriers using photo- 
CELIV are more reliable when Aj/j(0) -C 1, i.e. when r CT ^> t tT where t tr = d-Jlj [iV is the 
charge transit time through the layer. Under this approximation, the current density (CQ) 
becomes 

for 1(f) < d, where 

Ut) = \6d + 6de K \n* - 1) + U'tK 2 e K * («* - 3)1 , (10) 

and K*(t) = KyJU't/d. Unfortunately, equation does not provide any closed analytic 
expression for t max and Aj/j(0), but can be evaluated numerically. Figure M compares the 
apparent charge mobility calculated from such data using equation to the actual mobility 
at E = E* = U't max /d. These results were calculated for /i = 10~ 6 cm 2 /Vs, d = 100 nm 
and e r = 3, but are considered to be fairly general since they are independent of n at suffi- 
ciently low densities and invariant under the mutual transformation fiQ — > a/xo, U' —>■ all' 
for arbitrary a. It is obvious that significant errors in the apparent mobility occur at large 
k and large U'. We propose a simple improvement of the CELIV analysis by attributing 
the apparent mobility values to the extraction field redefined as E* = 0.65 ■ U't max /d. The 
corresponding relative error of the apparent mobility is also shown in figure M and stays 
within 20% in the relevant parameter regime. In order to test this approach, we numerically 
simulated CELIV transients in the Aj/j(0) 1 regime for various values of U', assuming 
d = 100 nm, e r = 3, /io = 10~ 6 cm 2 /Vs and k = 10~ 3 (m/V) 1 / 2 . Figure [9] compares the 
apparent mobility values determined according to equation (j3J) and associated with either 
choice of E*. We found that using E* = 0.65 • U't max /d generally gives better results for 
this type of field dependence. The experimental error can be further minimized by using 
an iterative procedure if the field-dependence of the apparent mobility indeed follows the 
Poole- Frenkel behaviour: (1.) measure CELIV transients at different U' to obtain a range 
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FIG. 8: Ratio of the apparent mobility yu app . and the actual mobility at the electric field of E* = 
U't max /d (left) and E* = 0.65 • U't ma ^/d (right) as function of the Poole- Frenkel parameter k. 
The apparent mobility was calculated from the current transients given by equation Q using 
Hq = 10~ 6 cm 2 /Vs, d = 100 nm and e r = 3 by applying equation ([3]) and is shown parametric in 
the voltage slope U' . 
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FIG. 9: Apparent charge mobilities calculated from simulated CELIV transients using equation ([3]) 
for field dependent charge mobilities as indicated by the solid line. The extraction field associated 
to determined mobility values was calculated as E* = U't max /d (squares), E* = 0.65 • U't max /d 
(circles) or using the iterative procedure as described in the text (triangles). 

of h and t max values via equation fl3]), (2.) determine preliminary parameters and 
from the measurement using E* = U't max /d, (3.) for each IT, calculate the theoretical 
CELIV transient from equation Qj, numerically evaluate (Aj/j(0)) th , t^ ax and /i th using 
equation and calculate 5 = dln(// th / n °^) /U't^^K^, (4.) associate each measured mobil- 
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ity value to the extraction field E* = 5U't max / d and determine optimized /ig and «W from 
this data, (5.) iterate the procedure by repeating steps 3 and 4 until the determined jj^ 
and k^ 1 ' stabilize. This procedure is only moderately complex but significantly enhances 
the accuracy of charge mobility determination, at least when Poole- Frenkel field dependence 
prevails. Figure [9] shows that the results of this iteration procedure accurately track the true 
field- dependent mobility used for the simulated CELIV experiments. 



V. CONCLUSION 



In this publication, we pointed out several difficulties that arise when applying the CELIV 
technique under the non-idealized conditions usually encountered in experiments. We based 
our analysis upon a rederivation of the original CELIV analysis, correcting for inaccuracies in 
the original publications that result in erroneous charge mobilities under conditions of high 
charge density. In the case of photogenerated charge carriers, we found a significant depar- 
ture from the equilibrium assumption of the original derivations. The impact of bimolecular 
charge recombination on CELIV transients and their analysis in this situation has not been 
considered up to now. We showed that high charge densities as typically used in the experi- 
ments leads to an artifical time dependence of the determined mobility values. We were able 
to relate the experimental results for a M3EH-PPV/CN-ether-PPV blend solar cell solely 
to this effect, showing that the charge mobility is actually constant shortly after photoex- 
citation. Sufficient information was provided to facilitate an interpretation of experiments 
under conditions of non-equilibrium charge carrier extraction. As another strong deviation 
from idealized conditions we investigated the effect of field-dependent charge carrier mobil- 
ities. We showed that association of the CELIV mobilities with the electric field present at 
the time of extraction current maximum leads to significant errors in the determined field 
dependence. An optimized choice of the correlated extraction field was introduced, which 
yields much lower errors compared to the standard approach. Additionally, we showed that 
under the presumption of a Poole-Frenkel type field dependence, an iterative procedure can 
be applied to determine the true mobility-field dependence in a precise way. 
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